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Abstract 

We present series study of using the method of super-symmetric quantum 
mechanics(S'C/ S'yQM) solving the spin-weighted spheroidal wave equation. In 
this paper, we obtain the first four terms of super-potential of the spin-weighted 
spheroidal wave equation in the case of s=l. These results may help summary 
the general form for the n-th term of the super-potentialwhich is proved cor- 
rect by means of induction. We finally compute the ground eigenvalues and 
ground eigenfunction. All the results may be of significative for studies of elec- 
tromagnetic radiation processes near rotating black holes and compute radiation 
reaction in curved space-time. 

Keywords: spheroidal wave equation, super symmetric quantum mechanics, super- 
potential, the ground eigenvalue and eigenfunction 
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1 Introduction 

The spin-weighted spheroidal functions are first defined by Teukolsky in the study of 
the perturbations of the Kerr black hole, which are indispensable in study of quantum 
field theory in curved space-time and black hole perturbation theory, etc. Its differential 
equations aret^l"^^' 
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where E is an eigenvalue to be found by imposing appropriate boundary conditions 
and the parameter s, which is the spin- weight of the perturbation fields could be 
s = 0, ±|, ±1, ±2 correspond to the scalar, neutrino, electromagnetic and gravitational 
perturbations respectively. And the boundary conditions requires is finite at ^ = 
0, TT. In series papers, the perturbation methods in the super-symmetric quantum 
mechanics(S'f/S'FQM) are used to study the spheroidal equations and many useful 
results are obtained. ^^'"^^^1 Here we continue to study the electromagnetic one, which 
is the spin- weighted equations in the case of s = 1. 

As done before, defining 0(^^) = and introducing an artificial small parameter /3 

au then Eq.(l) leadsW-t^^l 
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+ s + cos^9-2s(3 cosO- 
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+ E 



Thus, the corresponding boundary conditions being \l/|e=o 
in Eq.® is 



^1 



^ = (2) 
0. The potential 



V{9,f3,s) = - 



+ s + l3^ cos^ 9 - 2s(3 cos 6 



[m 



s cosy — i 



sin' 6 



(3) 



Super potential W is an important concept in SUSYQ, and it can be determined by 
potential with the relation as'^^' 

- W = V{e,(3,s) - Eo. (4) 
And the ground state wave function \E'o can be obtained from the super-potential by 



^0 



A^exp 



Wdz 



(5) 



Both the Eq.dl]) and Eq.(|2]) are both difficult to solve so that we look for the approxi- 
mation methods. In the paper, we will resolve the spin-weighted spheroidal equations 
by expanding the super-potential as a Taylor series in powers of /3 and the ground 
eigenvalue Eq can be also expanded by /jW"!^^] 



n=0 

oo 



(6) 
(7) 



n=0 



where n refers to the nth-order item of the series expansion, and m in EQ^nym, is referred 
to the parameter m in Eq.([2]) and the index means belonging to the ground state 
energy. Substituting Eqs.(l6])and ([7]) into Eq.(j4]), one could easily obtain 



1 [m + s cos 6 

Eo,0;m + -5 + r^—- 

4 sm^6' 



= MO) 
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W[-2WoWi = Eo,i;„,-2scos^ = /i(0) 

n-1 

W'^-2WoWn = E^,n,m + Y,^'^^n-k = > 3) 



k=l 



When s =1 , it is straight forward to obtain 



12 



0,0;m 



m +171 — 2 

1 + (m + i) cos 6* 
sinO 



From Eq.® it is easy to see that 

Wn{e) = {ts.n^-)-^sin~^'^^^eAn{e) = {l + cosefsin-^'^-''eA^{e) 

where, 

We introduce the following formulae to simplify calculation in Eqs. (fTO|) -( nA|) ^^^^ 



P(2m+ 1, 
where, 



cos 6 



2m + 2 



.k=0 



I{2m + l,k) 



(2m+2)(2m)---(2m+2-2fc) 'f > n 
(2m+l)(2m-l)-(2m-2fc+l) H ^ ^ U 

if A; < 



Here we give some recurrence relations between P[2m + 2n + 1, 6*) and P{2m 
for later use. Assuming 



P(2m + 2n + l,9) = Pi{2m + 2n + 1, ^) + P2(2m + 2n + 1, ^), 



where, 



Pi(2m + 2n + 1,1 
P2(2m + 2n+ 1,1 



— cos 9 
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— cos 6 
2m + 2n + 2 

— cos 6 
2m + 2n + 2 



sm 



^/(2m + 2n + 1, A;) 

fc=0 
m+n 

^ /(2m + 2n + 1,A;) 



2m+2n-2fc/: 



sm 



2m+2n-2k, 



.k=n+l 
m—1 



^/(2m + 2n + l,p + n + 1] 



sm 



2m- 2p 



p=0 



(2m - l)/(2m + 2n + 1, n + 1) 
2m + 2n + 2 



P(2m- 1, 



So Pi (2m + 2n + 1,9) contains the item like sin^™"'"^-' 9 with j = n — k > 0, which is 
the convergent item at 0, vr and P2{2m + 2n + 1, 9) can be expressed by P{2m — 1,9). 
Hence P{2m + 2n + 1, 9) can be rewritten as 

P(2m + 2n + l,9) = cost/ j(2m + 2n + 1, fc) 

(2.n-l)/(L°+2„+l.„+l) 

2m + 2n + 2 ^ ' ^ ^ ^ 

By using Eq.f ll7p . the calculation of Ai could be further simplified to be 



Ai{9) = j{Eo^i,^-2cos9){l-cos9)^sin^"'-^ 



9d9 



(2^0,1;^ + 4)P(2m - 1, - (Eo,i;m + 4)P(2m + 1, 
-sm 9 + 



m m + 1 

2m(Eo,i;m + 4)\ 
2^0,1;^ + 4 ^-^^^ J * P(2m - 1, 9) 

E0,l;m 'r2 . S%T?^^H Ep^,,^ + 4 . 

-sm 9 ^ 1 — cos 9stn 9 



m m + 1 2m + 1 

From Eqs.([IOD,(Il2]) the term P(2m - 1,6*) in ^1(6*) would make Wi{9) blow up at 
the boundaries ^ = 0, vr, that is to say unless the coefficient of the divergent term 
P(2m — 1, ^) in is zero, the eigen-function '^{9) could not finite at the boundaries. 
So 2Eo,i;m + 4 - = 0, that is 

^0,l;m = ^ (19) 

m + 1 

Equ.(19) can simplify Ai in the concise form 

Ai{9) = lJ^l-^[sin2^-2 + 2cos^] (20) 

Owing to Eq.f jTOj) . it is easy to obtain the Wi{9) 

Wi{9) = Ai{9){l + cos9)\sin9)-^"'-^ = ^sin^ (21) 

m + 1 

By the same calculation process, one could obtain W2{9) and 1^3(6') as 

, m(m + 2) cos 6* sin 6' m(m + 2) sin 6' 

^ " (2m + 3)(m + l)2 ~ (2m + 3)(m + 1)3 ^ ^ 

_ 2m cos g sing 2m sing msin=^g 

' ~ (m + l)4(2m + 3) ^ (m + l)5(2m + 3) (m + l)3(2m + 3) ^ '* 
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and the corresponding ii^o,2;m and £^o,3;m are 

+ 7m^ + 11m + 3 



Eo,2;m 



(m + l)3(2m + 3) 



^'^'■'"^ - ~(m + l)5(2m + 3) ^^^^ 

2 The general formula of n-th Wn for the super- 
potential 

In this section, we summarize the general formula for Wn from the results of Wi to 
as following 



[^1 



Wn{e) = COS ^ ^ an,k sin^^"^ ^ + 5Z ^"-'^ ^^^^''"^ ^ (25) 

A;=l fe=l 

where the coefficient an,k and fer^^fc need to satisfy 



an,k = if A; < 1 or /c > [|] 
bn',k = if A; < 1 or > [^] 



(26) 



Eq.( l25|l is a form different from the super-potential Wn in Ref.(9), that is 



n-l 



W„(^) = sin ^ ^ an,k cos^^+i n = 1, 2, • • • (27) 

k=0 

The next work is to prove its correctness by mathematical induction. Obviously 
Wn for = 1, 2, 3 satisfies Eq.( l25l) . Assuming that when N < n, Wn meet the general 
form, we need to verify Wn{0) also satisfy the general formula. Back to Eq.(l8]), taking 
Wkie) and Wn-kie) into ZZl W^W^-k, that is 

n-l [51+1 [f]+l 

fe=l p=2 p=2 



where hn^p and are constant coefficients: 



fl2l 



n— 1 p— 1 

hn,p ^ ^ ^ ^j ^^k,p—j^n—k,j ^k,p—l—j(^n—k,j ~l~ bk,p—jbn—k,j) 

k=l j=l 
n-l p-1 

S'njp ^ ^ ^ ^j (,(^k,p—jbn—k,j ~l~ bk^p—jdn—kj) (^9) 

fc=l j=l 
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In order to simplify the calculation process, the maximum limit of p can be written 
uniformly as [|] + 1, because of Eg. (1261) . we can get 



hn,p = 0, if p < 2 or p > [f ] + 1 
gn,p = 0, if p < 2 or p > [^^] 

Taking f„{9) = Eo^n;m + YlkZl WkWn-k into Eq.(IlI]) An{6) can be written as 



(30) 



0,n;m 



[fl+l [fl + l 

-E'o,n;m + ^ ^n,p siu^*'"^ 6* + COS 6* ^ ^(^^p siu^^"^ 
p=2 p=2 

sin^™ 6) 

2P(2m - 1, ^) - P{2m + 1,6)- 



(1 - coseysin'"'-^ede 



m 



gn,p sin2"+2*' ^{9n,p- hn,p) sin2'"+2p-2 ^ 



2m + 2p 



m + p — 1 



+ E 

p=2 

+ ^ [2{hn,p - 9n,p)Pi2m + 2p-3,e) + {2gn,p - K,p)P{2m + 2p-l,e)] (31) 

p=2 

We simplify the above result according to Eq. (fT7|) . that is 

An{e) 

= P{2m-l,e) xbi + Eo, 

[fl+i 



sin^™^ sin^"^^^ 

h cosfc*- - 

2m + 1 



m 



+ E 



p=2 
^+1 



gn,p Sin2-+2P ^ (^,p - /i„,p) sin2™+2p-2 ^ 



2m + 2p 

n,pj 

p=2 

M+i 

~l~ ^ ^ {2gri,p hn^p) 
p=2 



COS 



m + p — 1 



p~2 



2m + 2p-2 



— ^/(2m + 2p-3,p-2-j) 



2m+2jQ 



j=0 



— COS 6' 
2m + 2p 



p-i 



^/(2m + 2p- 1 ~j)sin^"'+'^^e 



3=0 



(32) 



where 6i is the coefficient of the term P(2m — 1 

[fl+i 

^1 = — ^ — ^ — + (2m - 1) 2^ 



2(m + l)i?o,n;m , ^ ^n,p " fl'n,p r/^ , o o i \ 

1/1™ i\ X .-r "^i (2m + 2p — 3,p — 1) 



2m + 



p=2 



m + p — 1 



+ (2m-l) 5^ 

p=2 



51+^ r 

2gn,p '''n,p 

2m + 2p 



/(2m + 2p - l,p) 



(33) 
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As similar analysis of ^1(6'), Wn{0) blows up at the boundaries 6 = 0, tt when bi ^ 
we must set the coefficient 61 of the term P{2m — 1,6) in An zero in order that the 
eigen-f unction \l/o(6') could finite at the boundaries. Thus 61 = gives 



E, 



0,n;Tn 



(2m + l)(2m- 1) 
2(m + 1) 



r]+l 



J2 ^"'P /"f /(2m + 2p - 3,p - 1) 



p=2 



m + p — 1 



(2m + l)(2m - 1) |- 2gn,p - K^ p 

2m + 2p 



2(m + 1) 



p=2 



-I{2m + 2p-l,p) (34) 



By the result of the Eg. (134!) . An{9) can be simplified as An{9) = RniO) + cos T„(6'), 
where 



Rn{e) 



with R 



n,0 



and 



^sin^- 
m 



p=2 



(^?n,p - K,p) sin2-+2p-2 ^ ^^^^ gi^2n^+2p ^ 



m + p — 1 



2m + 2p 



[fl+i 

J] i?„,p sin2'"+2p^ 

p=0 

and - ^"■^"^"■^ 



(35) 



^^-^ , and when p > 2, is 

2f7n,p+l 2/ln,p+l S'wjP 

2m + 2p 



(36) 
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m + p — 1 
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with 



-n,0 
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p=l 
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and when j > 1, Tnj being abbreviated as 



T 



p=j+2 



5'n,p iT"n,p 



m + p — 1 



I{2m + 2p-3,p-2- j) 



1+1 
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/(2m + 2p — l,p — 1 — j) 
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Due to Eq. fl2^ . one can get 



Rn,j = 0, ifj<-lorj>[^] 
T„,, = 0, ifj<-lorj>[f] 

According to Eq. (fTO|) . 

Wn{e) = [Rn{e) + cosOTnie)] (l + cos^)2sin-2'"-3^ = X„ + cos^y„ 

with 

X„ = [(2i?„ + 2T„,) - (i?„ + 2T„)sm2^] sin-2'^-=^^ 

= ^ [(2/?„j + 2Tn,j)sin^'~''9 - (i?„j + 2T„j)sm'^-^^ 

j=0 

[fl+i 
i=-i 

where 



2-Rn,j+l ~l~ '^Tjijj^i Rn,i 27~^j- 

2 

E 

p=j+i 



[(2m + 2p)/i„,p - 25(„J/(2m + 2p - 1 - j) gf, 



(38) 



(39) 



(40) 



(41) 



(42) 



m + j + l)(2m + 2p)(2m + 2p- 2) 2m + 2j 



(43) 
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And Yji can be expressed as 

y„ = [(2i?„ + 2r„)-r„sm2^]sin-2— 3^ 



[fl+i 

= J2ynjsin'^-'0, (44) 
j=-i 

from Eq.(35) and Eq.(38), we can get 

^ _ v^' (m + M2,n + 2j;)//„,, - 2(y„ J/(2//; + 'Ip - 1. p - 1 - j) 

(m + j + l)(2m + 2p)(2m + 2p-2) ^ ^ 

Obviously, in the case of j = —1,0, the items in X„ j and l^j of An{6) may cause 
Wn{9) divergent, then we need to checking whether these items are exist 

^n-i — 2i?„ + 2r„ = 

Xn,o — 2i?„_i + 2T„_i — Rnfi — 2T„_o = 

Yn,-1 — 2i?n,o + 2T„_o = 

Yn,o = 2i?„,i + 27;,i-7;,o = (46) 

Hence, all the items which can make Wn{0) divergent are not exist, so that Wn could 
be convergent in the boundary region. Wn{0) now becomes 

[fl+i [fl+i 
Wn{e) = J2 ^n,jSiri^^-'^9 + cose J^'^njsin^^-^e 

j=l 3=1 

= J]] fo„,iSm2^"^^ + cos^ ^ a„,/sm2'-^^ (47) 



1=1 1=1 



where 



^+1 



_ y _ _ (m + j)[(2m + 2p)hn.p - 2gn,p]I{2m + 2p - l,p - 1 - j) 
an,i - n,j- 2-^^ (m + j + l)(2m + 2p)(2m + 2p-2) 



2 



, ^ X = V [(2^ + 2p)/in,p - 2gn,p\I{2m + 2p - l,p - 1 - j) t/^j ^ 

(m + j + 1) (2m + 2p) (2m + 2p - 2) 2m + 2j ' 
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Note that the maximum hmit of / in Eq. ( H7|) is differ with the general formula Eq. ( 125|) . 
By the relation Eq. (HO|) we could analysis the maximum limit of / with the parity of n. 
When n is even, / gets the maximum / = [^] + 1 = ^ + 1, then 



On.f +1 = 2i?„^|+2 + 2T„^^+2 — ^n,f +1 = 

^n,^+l = 2i?„ "+2 + 2T„ "+2 — Rn,^ + l " 2T„ "+i = 



(49) 



but when ^ = [f ] = f , On.f and hn^Ti are both nonzero. Hence under the condition n is 



even Wn{0) becomes 



Wn{e) = cos 9 ^an,isin^^-^ 9 + J2bn,isin^^~^ 9, 



(50) 



When n is odd, / gets the maximum / = [^] + 1 



1=1 

then 



a„ n+l = 2i?„ n + 3 + 2T„ 11 + 3 — T„ n + 1 = 

2 ' 2 ' 2 ' 2 

&„,Il±i = 2i?„ n±3 + 2T„ n±3 - /?„^2i±l " 2T„ n±l 7^ 



(51) 



and when ^ = [f ] = o„ and b^ n^ are nonzero. While when n is odd Wn{9) 
has the new form 



Wni9) = cos^ J]a„,zsm2'^^0+ ^ 6„,zsm2'-i^, 



(52) 



z=i 



1=1 



with the help [f ] + 1 = [^^]- Through comparing the Eq.f lSOj) and Eq.( l52|) the super- 
potential Wn{9) can be rewritten as 



Wn{9) = cos9j2(^n,isin^^'^9+ J2^n,isin^^'^9 



(53) 



1=1 



This completes the suppose of the general formula Wn- That is to say, for any n > 1, 
Wn satisfy the form of Eq.(!25l). 

3 The ground eigenfunction 

From the super-potential 



W 



n=l 



cos 9 ^ an,k sin^'' ^ ^ + X] ^"'^ ^^^^^ ^ ^ 



k=l 



k=l 



(54) 
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the ground eigenfunction becomes 
\i/Q = N exp 

= iVexp 



- / Wdz 

1 + (m + i) cos6' 



sinO 



N{l-cose) sin^+s^exp 



n=l 



(i^ - ^ /3" / (cos e J2 an,k + Yl ^^'k) sin^'^-i ede 



k=l 

sin^^ 6 



k=l 



_ rii+ii 



n=l \ fc=l 



k=l 



(55) 



and the ground eigenvalue is 

oo oo 



(56) 



n=l 



n=l 



Eo^n;m has bccu obtained by Eq. (!Ml) . 

References 
References 

[1] Teukolsky S.A., Phys. Rev. Lett. 29, 1114 (1972); Astrophys.J. 185, 635 (1973). 

[2] Berti E, Cardoso V and Casals M 2006 Phys. Rev. D73 024013 

[3] Marc Casals and Adrain C.ottewill, High frequency asymptotics for the spin- 
weighted spheroidal equation. Phys. Rev. D , 2005,71:064025; 

[4] Tian Guihua, Shuquan Zhong, Solve spheroidal wave functions by SUSY method, 
arXiv: 0906.4685v2 [gr-qc] 30 Jun 2009. 

[5] Tian Guihua, Shuquan Zhong, New investigation for the spheroidal wave functions, 
arXiv: 0906.4687v2 [gr-qc] 30 Jun 2009. 

[6] Tian G H, 2005 Chm. Phys. Lett. 22 3013 

[7] Tian Guihua and Zhong Shuquan 2010 Chm. Phys. Lett. 27 040305. 

[8] Tian G H and Zhong S Q 2009 Arxiv: 0906.4687 V3: Investigation of the recur- 
rence relations for the spheroidal wave functions , preprint. 

[9] Tang W L and Tian G H 2010 Solving the spheroidal wave equation with small c 
by SUSYQM method, accept by Chin. Phys. B, preprint 



11 



[10] Zhou J, Tian G H and Tang W L 2010 The spin-weighted spheroidal wave functions 
in the case of s — ^, J. Math. Phys. (submitted) 

[11] Li K, Sun Y, Tian G H and Tang W L 2010 A tentative study of the spheroidal 
wave function in the case of s=2 , accept by Chin. Sci. G, preprint (In Chinese) 

[12] Dong K Tian G H and Sun Yue 2010 Study of the Spin-weighted Spheroidal Wave 
Equation in the Case of s = 1/2 Chin. Phys. Lett, (submitted) 

[13] Tian G H 2010 Chin. Phys. Lett. 27.030308 

[14] Cooper F, Khare A and Sukhatme U 1995 Supersymmetry and Quantum Mechan- 
ics 

[15] Gradsbteyn I.S., Ryzbik L.M. 2000 Table of integrals, series, and products. 6th 
edition, Elsevier ( Singapore) pte. Ltd. 



12 



